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ABSTRACT 

We investigate the spiral galaxy NGC 5394, which is strongly interacting with the 
larger spiral NGC 5395 (the pair is Arp 84), using optical integral-held spectroscopy 
from the CALIFA (Calar Alto Legacy Integral Field Area) survey. Spatially-resolved 
equivalent widths, emission-line ratios and kinematics reveal many features related 
to the interaction, which has reshaped the galaxy. Ha maps (with other diagnostic 
emission lines) show a concentrated central (r < 1 kpc) starburst and three less 
luminous star-forming regions (one knot far out in the northern arm), and we estimate 
the dust-corrected total star-formation rate as 3.39 Mgyr -1 . However, much of the 
galaxy, especially the outer tidal arms, has a post-starburst spectrum, evidence of a 
more extensive episode of star-formation a few x 10 8 yr ago, triggered by the previous 
perigalacticon. The [NII]6584/Ha ratio is high in the nucleus, reaching 0.63 at the 
centre, which we interpret as related to high electron density ( n e ~ 750 cm~ 3 from 
the [SII] |||y ratio). We find a central region of strong and blueshifted Nal(5890,5896) 
absorption, indicative of a starburst-driven outflow from the nucleus at an estimated 
velocity ~ 223 km s -1 . The CALIFA data also show an annular region at radii 2.25-4 
kpc from the nucleus, with elevated ratios of [Nil], [OI]6300 etc. to the Balmer lines - 
this is evidence of shock excitation, which might be the result of interaction-triggered 
gas inflow. 

Key words: galaxies: individual: NGC 5394, galaxies:interactions, galaxies: evolu¬ 
tion, galaxies: starburst 


1 INTRODUCTION 

Mergers of two spiral galaxies are one of the most impor¬ 
tant processes in galaxy evolution, and typically occur in a 
sequence of stages over several hundred Myr, with successive 
close passages (perigalactica) of the pair, multiple bursts of 
star-formation, and complex morphological transformations 
preceding the final coalescence into a single galaxy (which 
could become an elliptical or lenticular). Interacting galax¬ 
ies at any stage can be well studied by integral field spec¬ 
troscopy (IFS), which may reveal both the current activity 
and recent history within each resolution element. 

NGC 5394 and NGC 5395 are an interacting close pair 


* nathan.roche@astro.up.pt 


of spirals, which appear to be destined to merge. This well- 
known system was included in the Arp (1966) catalog of pe¬ 
culiar galaxies under the single name of Arp 84, and studied 
further by Arp (1969). On account of its shape it has also 
been called the Heron Galaxy, with the larger NGC 5395 be¬ 
ing the ‘body and wings’ and the two-armed, northern NGC 
5394 the ‘neck, head and beak’; this is well-depicted on the 
Gran Telescopio Canarias gallery image (Fig 1). 

NGC 5394, a barred spiral seen almost face-on and 
classed as SB(s)b-pec, appears relatively symmetric, but its 
structure is unusual and greatly modified by the interaction. 
Within its compact disk it has a small, bright nucleus and 3 
arms (double on the west) and beyond the disk (radii r > 13 
arcsec out to ~ 40 arcsec) a distinct and much larger system 
of two long arms extending north and south, the southern 
arm reaching NGC 5395. 
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The Arp 84 system is described in some detail by Kauf¬ 
man et al. (1999), hereafter K99, who observed in Ha with a 
Fabry-Perot. NGC 5394 and 5395 have a mass and luminos¬ 
ity ratio about 1:4, projected separation 28 kpc, and rotate 
anticlockwise and clockwise respectively, so that their arms 
are trailing. The interaction (an orbit approximately in the 
sky plane) is prograde with respect to the rotation of NGC 
5394, but retrograde and inclined for NGC 5395. NGC 5394 
has strong Ha emission (indicating star-formation) from the 
nucleus (2.2 x 10 -13 erg cm -2 s -1 from Keel et al. 1985) and 
more from a region in the SW of the disk (the inner western 
arm) and two bright spots on the inner and outer north¬ 
ern arm. K99’s radio observations give the neutral hydrogen 
(HI) mass in NGC 5394 as 7.3 x 10 s Mg, and the CO data 
(tracing molecular hydrogen) show a dense concentration 
(> 10 9 Mg) of Ht in the nucleus. K99 describe NGC 5394 
as ‘post-ocular’ in that it resembles the ‘eye-shaped’ inter¬ 
acting galaxies IC 2163 and NGC 2535 but is at a slightly 
later stage where ‘the eye-shaped oval has evolved into inner 
spiral arms and most of the gas has fallen into the central 
region, where it fuels a starburst’. K99 reported some evi¬ 
dence for a low-velocity outflow from asymmetry of the Ha 
line (a possible blueshifted component) and small (~ 10 km 
s -1 ) distortions in the Ha velocity map (NE of the nucleus); 
unconfirmed but will be investigated here. 

Kaufman et al. (2002), hereafter K02, performed further 
CO observations and comparison with several other pass- 
bands. They estimate NGC 5394 is inclined 15° ± 2 to the 
sky plane, discuss evidence (e.g. 60jrm flux) the Ha flux is 
reduced by as much as an order of magnitude by internal 
dust extinction, and on this basis estimate the total (dust- 
corrected) star-formation rate (SFR) as 6±2 Mgyr -1 (with 
no evidence of an active galactic nuclei - AGN - contribu¬ 
tion). The star-formation closely traces the HI and H 2 dis¬ 
tributions - all three are concentrated in the central kpc, 
but also the disk is ‘lopsided, with more HI, CO, and Ha 
emission coming from the western or southwestern side’, and 
there was speculation that this resulted from a west-vs-east 
time lag in the formation of disk arms and subsequent in¬ 
fall of gas to the nucleus. They estimate a molecular/atomic 
ratio M(H 2 )/M(HI) = 2.5-2.7 (even with a lower CO-H 2 
conversion factor applied for starbursts), very high com¬ 
pared to 0.86 for NGC 5395, which contains much more HI 
(> 10 10 Mq), and the average ratio of 0.28 for spiral galaxies 
(Casoli et al. 1998). 

Puerari, Valdez-Gutierrez & Hernandez-Lopez (2005) 
studied the morphology of both galaxies in the near-infrared 
(JHK), and fit (in the A-band) a bulge profile with r e // = 
3.12 arcsec plus a disk of scale length h = 6.35 arcsec for 
NGC 5394, which is ‘more compact than normal galaxies’. 
From the disturbed, Cartwheel-like structure of NGC 5395 
they conclude the system is undergoing a strong interaction, 
a ‘collision rather than a grazing encounter’. Kaneko et al. 
(2013) again study the system in CO and depict the contrast 
between the compact, nucleated, H 2 -rich NGC 5394 and the 
ring-like and Hi-rich 5395. 

Smith et al. (2007) observed the system from 3.6 to 
24/nn in a Spitzer survey of Arp-catalog galaxies: most emis¬ 
sion from NGC 5394 is from the compact nucleus and the 
Spitzer images also show the prominent bright spot (Ha 
source) out in the northern arm. Higdon et al. (2014) de¬ 
scribe this hotspot (‘Nl’) as an ‘intergalactic star-forming 



Figure 1. Gran Telescopio Canarias public gallery im¬ 
age of Arp 84, NGC 5394 (north) 5395 (south). From 
www.gtc.iac.es/multimedia/media/GTC_ARP84.jpg 

object’ with a 8/mi flux of 2.626 ± 0.015 mjy, ~ 1% that 
of the whole galaxy (279.8 ± 3.8 mjy), and in the mid-IR 
it is also detected in [Nell], [Neill], [Sill] and polycyclic 
aromatic hydrocarbons (PAHs). 

Note that the SFRs estimated from observable quanti¬ 
ties will generally depend on the assumed initial mass func¬ 
tion (IMF), e.g. the SFR from K02 assumed a Salpeter IMF, 
whereas equivalent SFR estimates based on the Chabrier or 
Kroupa IMFs, with fewer low-mass stars, will be 40-45% 
lower. Lehmer et al. (2010) estimated a (Kroupa IMF) SFR 
of 9.9 Mgyr -1 for the whole Arp 84 system, simply from the 
infra-red luminosity of 10 11 Mgyr -1 , and found its X-ray 
luminosity from Chandra data consistent with the typical 
L x — SFR relation of luminous infra-red galaxies (LIRGs), 
meaning it can be accounted for by X-ray binaries without 
requiring AGN. Howell et al. (2010) estimated a higher SFR 
of 20.99 Mgyr -1 by including the unobscured component 
traced by far-UV luminosity ( GALEX) as well as the infra¬ 
red luminosity and modelling with the Salpeter IMF. 

Lanz et al. (2013) combined fluxes in many bands from 
FUV [GALEX) to FIR ( Herschel ) for a sample of interact¬ 
ing galaxies and fitted the SEDs with a program MAGPHYS, 
which assumed a Chabrier IMF and gives a SFR averaged 
over the past 100 Myr. For NGC 5394 and 5395, they esti¬ 
mate SFRs of 1.54 and 3.69 Mgyr -1 , and found some differ¬ 
ence in their SEDs, with NGC 5394 having relatively more 
12-60/im flux and higher dust temperatures (23.2K in the 
cold component) in comparison to 5395 (with 19.7 K) and 
the non-interacting control sample (about 19K). This is con¬ 
sistent with NGC 5394 experiencing a strong (‘stage 4’), 
and late-stage interaction, with star-formation more cen¬ 
trally concentrated than in the more massive companion. 
Stellar masses are estimated as 4.79 and 15.14 x 10 lo Mg. 

In this paper we investigate the spatially-resolved prop- 
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erties and history of NGC 5394, especially the effects of 
the interaction, using new integral field spectroscopy from 
the Calar Alto Legacy Integral Field Area (CALIFA) survey 
(Sanchez et al. 2012; Husemann et al. 2013). CALIFA has al¬ 
ready been the basis for a number of papers on global prop¬ 
erties and on in-depth studies of particular galaxies, with 
which we can compare - e.g. the famous ‘playing Mice’ inter¬ 
acting pair NGC 4676A/B (Wild et al. 2014), and studies of 
metallicity gradients in typical spirals and mergers (Sanchez 
et al. 2014; Barrera-Ballesteros et al. 2015). We will also 
compare with models of emission-line ratios in star-forming 
galaxies (e.g. Kewley et al. 2013) and in shocks (Allen et al. 
2008); and simulations of galaxy mergers (e.g. Barnes 2004, 
Struck et al. 2005, Torrey et al. 2012). 

From the CALIFA catalog, the J2000 co-ordinates of 
NGC 5394 are R.A. 13:58:33.19 and Dec +37:27:13.11 and 
the heliocentric optical recession velocity is 3457 km (z = 
0.01153); the velocity relative to the local standard of rest 
(Virgo cluster) is given as a slightly greater 3688 km s -1 
(corresponding to z = 0.01230). For Ho = 70 km s _ 1 Mpc^ 1 , 
11 m = 0.27, 12a = 0.73 (assumed throughout) the Virgo 
velocity corresponds to a distance modulus 33.64 mag (the 
value listed in the CALIFA catalog) and an angular diameter 
distance 52 Mpc, so that l arcsec subtends 252 parsec. In 
Section 2 of this paper we describe the data, in Section 3 
the distribution of star-formation activity, in Section 4 and 
5 diagnostic line ratios and metallicity, in Section 6 the line- 
of-sight kinematics and and Section 7 discuss the history of 
the interaction (summarizing in Section 8 ). 


2 CALIFA OBSERVATIONS 

Our data set for NGC 5394 comes from the first CALIFA 
public data release (DRl; Husemann et al. 2013), which 
consists of optical datacubes (image x wavelength) for the 
first 100 (of what will be eventually 600+) galaxies in the 
survey, now updated with Data Release 2 (Garcfa-Benito 
et al. 2015). CALIFA observations are obtained with the 
integral-held spectrograph Potsdam MultiAperture Spec¬ 
trophotometer (PMAS)/PPAK mounted on the 3.5m tele¬ 
scope at the Calar Alto observatory. The PPAK unit con¬ 
tains 331 densely packed optical fibres arranged in a hexago¬ 
nal area of 74 x 65 arcsec, to sample an astronomical object 
at 2.7 arcsec per fibre, with a filling factor of 65% due to 
interfibre gaps. 

The CALIFA galaxy sample is selected from the Sloan 
Digital Sky Survey (SDSS) DR7 catalog, simply on the basis 
of redshift (0.005 < z < 0.03) and apparent size (angular 
isophotal diameter 45 < D 25 < 80 arcsec), so that each 
galaxy fitted within, while filling most of, a single held-of- 
view. This means the survey covers a very wide range of 
morphological type, absolute magnitude, stellar mass, star- 
formation rate, colour, etc. (Walcher et al. 2014). NGC 5395 
was not targeted in CALIFA as its apparent size is more than 
twice NGC 5394 (radii Rk 2 o = 34 and 71 arcsec; Kaneko et 
al. 2013) and exceeds a single held-of-view, so this analysis 
is limited to the more compact galaxy. 

Each galaxy was observed with two spectrograph config¬ 
urations, (i) the V500 grating covering the wavelength range 
3745-7500A with a spectral resolution of 6.0A (full-width at 
half-maximum - FWHM), and (ii) the V1200 grating cover¬ 


ing 3650-4840A with a spectral resolution of 2.3A (FWHM). 
Exposure times are 900 sec on each galaxy with V500 and 
3 x 600s with V1200, and the respective 3<r limits for de¬ 
tection of a narrow emission line are 1.0 and 0.6 x 10 1 ' 
ergs cm~ 2 s _1 . Initial processing of these data generated two 
(V500 and V1200) data cubes for each galaxy, which are 
supplied on the CALIFA website (califa.caha.es) calibrated 
in units of A and 10 -16 erg cm - 2 s - 1 A -1 and rebinned into 
a grid of l arcsec spatial pixels. As the spectrograph fibre di¬ 
ameter is 2.7 arcsec, the effective spatial resolution is lower, 
with FWHM from 3 to 4 arcsec. The cubes are 78 x 73 
(spatial) x 1877 (spectral) with dA = 2 . 0 A for V500, and 
78 x 73 x 1701 with dA = 0.7A for V 1200 . 

In this paper we concentrate mostly on the V500 spec¬ 
tra which contain the lines of interest. We have analyzed the 
datacubes using a locally developed software pipeline called 
PORT03D (Papaderos et al. 2013, Gomes et al. 2015), which 
effectively separates out the stellar and nebular (emission¬ 
line) components of each spatial pixel’s spectrum. The stel¬ 
lar component is fitted with a set of model stellar spectra, 
with a wide range of ages, using the STARLIGHT package 
(Cid Fernandes et al. 2013). The total stellar mass is esti¬ 
mated from the PORT03D fit as 4.6 x 10 10 Mg (for Salpeter 
IMF). For a list of spectral lines, the program generates 2D 
maps of emission-line flux, equivalent width, velocity offset 
(red or blueshift), and the equivalent widths for any stellar 
absorption lines (for Balmer lines H/Lyd there may be nar¬ 
row emission within broader absorption; as was seen for this 
galaxy by Arp 1969). Further maps show stellar properties 
estimated from the model fits, such as mean stellar age and 
metallicity (mass and luminosity weighted). One line map, 
for [OI]6300, appeared mildly contaminated by the sky line 
[OI]6364 (de-redshifted to 629lA), so to exclude this a new 
flux map was computed summing 6297-6307A only. 


3 MAPPING THE STAR-FORMATION 

Figure 2 shows the CALIFA image of NGC 5394 in the con¬ 
tinuum at 6390-6490A, just blueward of Ha, where it ap¬ 
pears near-symmetrical, but the emission lines will reveal 
a much more complex picture. The Ha emission line is the 
best tracer of the current star-formation rate (SFR). CAL¬ 
IFA maps of the Ha flux and equivalent width (EW, here in 
the rest-frame) show (Fig 3) the star-formation is concen¬ 
trated in four distinct regions, the nucleus (r < 4 arcsec or 
1 kpc), a region in the SW of the disk (identified by K99 
as the inner of two western spiral arms) and two compact 
‘hotspots’ in the inner and outer northern tidal arm (see 
also Fig 4 of K99). Peak Ha equivalent widths are 52.8, 70.3 
and 54.1 A for the nucleus, SW disk and inner-N hotspot, 
with the highest of all, 106.7A for the far-N hotspot, but 
the nucleus produces by far the greatest Ha luminosity. 

Arp (1969) described an ‘unusually sharply defined 
string of condensations coming out of the companion on 
the inside of the northern arm’. This ‘string’ is visible on 
the GTC image, and the inner-N and far-N Ha-emitting 
hotspots can be identified with the knots at either end. 
Probably the intermediate knots are recent sites of star- 
formation, which now continues only at the extremities. 
There is lower intensity Ha emission distributed over the 
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NGC 5394 Log-scale continuum 6390-6490 Ang. 
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Figure 2. NGC5394 red-band continuum map from CALIFA 
data. North is up and East is left, as in all figures. 


whole disk (out to r ~ 12 arcsec), but almost none from the 
tidal arms outside of the two hotspots. 

To illustrate some of the diversity within this galaxy, 
Fig 4 compares small-aperture CALIFA spectra of the nu¬ 
cleus, the far-N hotspot (averaged over the central 5 pixels 
of each) and a region in the southern arm (averaged over 
14 pixels). Wavelengths 5508-5522A, 6222-6236A and 6285- 
6295A (galaxy restframe) are severely affected by strong sky 
OI lines (observed at 5577.3, 6300.3 and 6363. 8A), and are 
masked out in the plots and model fitting. The nucleus has 
the strong emission lines expected for a starburst, although 
[OIII]5007 is relatively weak and the continuum relatively 
red (i.e. redder than flat in F\). The two outer regions are 
bluer and the N hotspot has strong emission lines, while the 
S arm is not only passive but may be post-starburst on the 
basis of strong Balmer absorption lines (EW > 5A). The 
nucleus has notable differences in emission-line ratios com¬ 
pared to the hotspot, with higher [NII]6584 /Hq and lower 
[SII] 6717/ [SII] 6731. 

SFRs locally and for the whole galaxy can be estimated 
by summing the Ha fluxes, but there is evidence (K02) 
the emission lines are substantially attenuated by internal 
dust. For star-forming HII regions the intrinsic flux ratio of 
Ha/H/3 (when corrected for stellar absorption as here) is ex¬ 
pected to be 2.86 (from ‘case B’ models with T e ~ 10 4 K), 
but the observed ratio is increased in proportion to the dust¬ 
reddening. Using the extinction law of Calzetti et al. (2000), 
supported by the more recent results of Wild et al. (2011), a 
dust-corrected Ha flux is estimated for each individual pixel 
as F(Ha) CO rr ~ F(Ha) oba [Ha/2.86H/?] 2 ' 615 . Here this dust 
correction is only applied for pixels where the statistical er¬ 
ror cr(Ha/H/3) < 1.0, which includes all of the disk and the 
star-forming regions in the northern arm (the other pixels 
are included in the Ha flux summations but uncorrected). 
Figure 5 shows a map of the Ha dust extinction, which is 
greatest for the nuclear region (2.32 mag). 

To depict the IFS data (without defeating the object 
by too much spatial binning) can be challenging and three 
methods are used here, 2D maps, radial profiles averaged 


NGC 5394 Ha Emission Equivalent Width 
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Figure 3. Map of (log scale) Her emission-line flux arcsec 2 
(above) and equivalent width in A (below). 


in circular annuli, and we also divide the galaxy into 7 re¬ 
gions, shown on Fig 6. Table 1 gives the observed and dust- 
corrected Ha fluxes, summed over each of the 7 regions and 
for the whole galaxy (for the integrated Ha the estimated 
extinction is 1.98 magnitudes). Fig 7 shows the integrated 
CALIFA spectra of each region (note they have different 
areas), which will be model-fitted in a later Section. 

From the dust-corrected fluxes, Ha luminosities in erg 
s -1 are calculated using the distance modulus 33.64 mag. 
A star-formation rate in solar masses (Mg) per year can 
then be estimated as 4.4 x 10 -42 L(Ha) (e.g. Sobral et 
al. 2014); this conversion factor is for the Chabrier IMF. 
We estimate 3.39 Mgyr -1 for the whole galaxy. For the 
Salpeter IMF (with more low mass stars) the SFR would 
be 7.9 x 10 -42 L(Ha) (Kennicutt 1998) giving 6.09 Mgyr -1 , 
in agreement with the K02 estimate 6 ± 2Mgyr -1 . In com¬ 
parison, for the Mice galaxies (each larger and about twice 
as massive as NGC 5394), Wild et al. (2014) estimated SFRs 
6.2 Mgyr -1 for NGC 4676A and 2.3 Mgyr -1 for NGC 4676B 
(correcting for dust using UV and mid-infrared fluxes). 
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Region 

F(Ha)obs 
erg cm -2 s —1 

F(Ho;)corr 
erg cm -2 s _1 

L(Ho;)corr 
erg s —1 

SFR corr 

M 0 yr _1 

[Nil] 

Ha 

[OIII] 

H/J 

[SHI 

Ha 

[OI] 

Ha 

[sn]§M 

Nucleus 

1.981 x 10~ 13 

1.681 x 10" 12 

5.749 x 10 41 

2.53 

0.592 

0.273 

0.288 

0.023 

1.069 

SW disk 

3.776 x 10~ 14 

1.547 x 10~ 13 

5.290 x 10 4 ° 

0.23 

0.452 

0.262 

0.296 

0.031 

1.375 

Outer disk 

9.414 x 1CT 14 

3.603 x 10" 13 

1.232 x 10 41 

0.54 

0.551 

0.407 

0.343 

0.043 

1.351 

Inner-N hotspot 

5.273 x lCT 15 

1.154 x 10" 14 

3.946 x 10 39 

0.017 

0.466 

0.425 

0.289 

0.033 

1.301 

Far-N hotspot 

7.257 x 10" 15 

2.031 x 10" 14 

6.946 x 10 39 

0.031 

0.408 

0.280 

0.270 

0.019 

1.526 

N arm 

1.161 x 10~ 14 

1.578 x 10" 14 

5.396 x 10 39 

0.024 






S arm 

8.861 x 1CT 15 

1.088 x 10" 14 

3.721 x 10 39 

0.016 






Whole Galaxy 

3.629 x 10~ 13 

2.242 x 10" 12 

7.710 x 10 41 

3.39 

0.562 

0.368 

0.312 

0.033 

1.191 


Table 1. Integrated Ha luminosities for each of 7 regions of the galaxy NGC 5394 (Fig 6); fluxes and luminosities corrected for dust 
using the Balmer decrement; star-formation rates derived from the corrected L(Ha) (Chabrier IMF); and several important line ratios 
from integrated, uncorrected fluxes (not shown for tidal arms as the line fluxes are so low). 


NCG 5394 (spectra of regions) 




Wavelength Angstroms 

Figure 4. CALIFA V500-grating spectra in 3 different regions of 
NGC 5394, the nucleus, the far-north star-forming hotspot (av¬ 
eraged over central 5 pixels), and the (post-starburst) southern 
tidal arm. 


We estimate the far-N hotspot accounts for 0.9% of 
the galaxy’s SFR, the same as its fraction of the 8/tm flux 
(2.626/279.8 mjy) in Spitzer observations (Higdon et al. 
2014). Such features are not unexpected e.g. the ‘ocular’ IC 
2163 has at least 9 star-forming clumps spaced around the 
‘eye’, with similar Ha and 8 pm fluxes (Kaufman et al. 2012). 
In NGC 5394, star-formation is now concentrated in the nu¬ 
cleus, with 75% at r < 1 kpc. Fig 8 shows the observed and 


Ha dust extinction (from Ha/HfS ratio) 
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Figure 5. Her dust extinction as estimated from the Balmer 
decrement Ha/H/3 using the Calzetti et al. (2000) law. 


NGC 5394 star-forming and other regions map 
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Figure 6. Map of regions the galaxy is divided into for the anal¬ 
ysis in Table 1: (1) S Arm (2) N arm (3) Outer disk (4) Inner-N 
hotspot (5) Far-N hotspot (6) SW disk (7) Nucleus 
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NCG 5394 (integrated spectra of regions) 



Figure 7. CALIFA V500-grating spectra integrated over 7 re¬ 
gions of NGC 5394 as shown in Fig 6, on a log scale. 


NGC 5394 H8 Absorption Equivalent Width 
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NGC5394: Ha Radial Profile 
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Figure 8. Radial profile of Ha emission, centred on the nucleus, 
as observed and as corrected for dust using the Ha/H (3 ratio. 
Also shown, at r < 10 arcsec, the profile of the hotspot in the 
northern arm, and for comparison, the shallower profile of the 
red continuum (6390-6490A) in flux A A 


dust-corrected Ha flux profile, averaged in circular annuli 
centred on the nucleus (the peak at r = 29 arcsec is the far- 
N hotspot), and also profiles centred on the far-N hotspot 
itself out to r = 10 arcsec only. The nucleus and hotspot 
are very different in luminosity but are concentrated star- 
formation regions and have similarly steep profiles. In the 
continuum, the profile of the galaxy is shallower, closer to 
an exponential, without the steepening at r < 10 arcsec. The 
central star-formation density is very high (higher than in 
NGC 4676A), with an integrated Ha flux 1.03 x 10~ 13 (ob¬ 
served) or 9.39 x 10 -13 (corrected) in the central 14 pixels 
(0.889 kpc 2 ); the latter corresponds to 1.6 Mgyr^kpc -2 . 


4 LINE RATIOS 

We examine a number of other diagnostic lines and ratios. 
First, the equivalent width of the H<5 stellar absorption line 
(Fig 9) is high (6-9A) over most of the galaxy area, espe¬ 
cially the tidal arms. An EW a bs(H<5) > 5A is considered 


Figure 9. Map of H(5 absorption-line equivalent width. 


diagnostic of a post-starburst region or galaxy (e.g. Goto 
2008) and is caused by an enhancement in the fraction of 
A-type stars, which on the basis of stellar lifetimes implies 
they formed in a period of increased star-formation between 
100 Myr and 1 Gyr ago. Only in the active star-forming 
(Ha-emitting) regions is EW(H<5) reduced to a more typical 
3-4A, presumably because the A-type spectrum is mixed 
with even younger OB stars. This shows NGC 5394 is very 
much a post-starburst galaxy with the earlier burst more 
spatially extended than the current activity. 

Of much importance in identifying AGN and other 
sources of high excitation, and estimating metallicity, are the 
line ratios [NII]6584/Ha (the N2 index) and [OIII]5007/H/1 
(the 03 index). These and the other line ratios considered 
here (except Ha/H/I) are almost insensitive to dust because 
the two wavelengths are very close. For the galaxy as a whole 
03 = 0.368 and N2 = 0.562. Fig 10 and 11 show the ratio 
maps, and Fig 12 the radial profiles (centred on the nucleus, 
and for comparison at the far-N hotspot). Many pixels, gen¬ 
erally in the outer galaxy where emission lines are weak or 
absent, give very noisy or meaningless ratios, and these maps 
and profiles include only pixels where the ratio measurement 
is > 3cr. 

The 03 ratio is low at 0.25-0.3 within the main star¬ 
forming regions, with no increase at the centre (only 0.254 
for the central 4 pixels), so gives no evidence for an AGN. In 
contrast, N2 is strongly peaked at the nucleus, 0.619 for the 
central 4 pixels and 0.629 at the highest pixel, but falls to 
0.52 at r = 6 arcsec, and 0.41 in the far-N hotspot, which is 
more typical for star-forming regions. However, at the outer 
edge of the disk the mean N2 increases again to > 0.6, and in 
many pixels here it is even higher than at the centre, with 
ratios 0.65-0.85. At the same radii there is an increase in 
03, and to some extent in [SII]/Ha and [OI]/Ha. We find 
an outer region of enhanced line ratios forming an irregular 
annulus, mostly between r ~ 9 and 16 arcsec. The N2 ratios 
in both the nucleus and this outer region exceed the log 
(N2)=-0.25 upper limit of typical star-formation (e.g. the 
Kauffmann et al. 2003 divide). 

CALIFA data for the Mice galaxies (Wild et al. 2014) 
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Figure 10. Map of [NII]/Ho ratio (N2 index). 


Figure 13. Map of [OI]6300/H« ratio. 
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Figure 11. Map of 0[III]5007/H/3 ratio (03 index). 


Figure 14. Map of the S[II]6717+6731/Ha ratio. 


NGC5394: Radial Profiles of N2 and 03 line ratios 
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Figure 12. Radial profiles of the line ratios N2 and 03 (centred 
on the nucleus), and for comparison the profiles centred on the 
far-N hotspot. 
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showed large biconical regions in both disks with N2 > 
—0.25 dex (0.56) up to N2 ~ 1, with spatially matching 
high values for the ratios [OI]6300/Ha (with 0.1-0.25) and 
[SII]/Ha (with 0.6-1.0), and (in 4676A) 03 (which ranged 
from 1-4). These high ratios and the bicone morphology 
could be explained by high-velocity (~ 350 km s” 1 ) shocks 
with ‘ionized precursor’, from a super wind, an outflow 
driven by a central starburst. If the nuclear starburst in 
NGC 5394 causes similar shock ionization, this would ap¬ 
pear directly in front of the nucleus due to the face-on ori¬ 
entation. However, at the nucleus, 03 and the other two key 
ratios [OI]6300/Ha and [SII](6717+6731)/Ha (Fig 12 to 15) 
do not show high values or any peak relative to larger radii, 
and all three (in contrast to the N2 ratio), remain within 
the range expected for HII regions, and so give no evidence 
of shocks (or AGN) here. 

The flux ratio of the two lines in the [SII] 6717,6731 
doublet is a useful diagnostic of electron density (n e ) in neb- 
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NGC5394: Radial Profiles of S2 and [Ol]6300/Ha line ratios 
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NGC5394: Radial Profiles of [SIIJ6717/[SII]6731 (electron density diagnostic) 
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Figure 15. Radial profiles of the line ratios OI/Ha and 
[SII](6717 + 6731)/Ha in the disk, and centred on the far-N 
hotspot. 


Figure 17. Radial profiles of the line ratios [SII]6717/[SII]6731 in 
the disk, and centred on the far-N hotspot. An n e scale is shown 
for T e = 10000if. 



Figure 16. Map of [SII]6717/[SII]6731 ratio, sensitive to electron 
density n e 
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Figure 18. Plot of the line ratios N2 vs. [SII]6717/[SII]6731 in 
the strongly star-forming regions, showing that they are corre¬ 
lated (the correlation coefficient is 0.8375) forming a sequence of 
increasing n e . An n e scale is shown for T e = !0000 A. 
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ulae, most sensitive at n e ~ 10 2 - 10 4 cm" 3 (with also a small 
dependence on temperature, see e.g. Canto et al. 1980). In 
NGC 5394, the [SIIJ6717/6731 ratio (Fig 16 and 17) shows 
a strong monotonic gradient at r < 6 arcsec, falling from 
1.4 in the outer disk (r ~ 6-12 arcsec) to 0.97 at the cen¬ 
tral 5 and 0.949 in the central one pixel. For the calibration 
in iraf.stsdas.nebula.temden, assuming T e = 10000K, this 
signifies a rise in electron density from n e — 27 cm -3 to a 
central 750 cm" 3 . The higher [SII]6717/6731 ratio of 1.55- 
1.6 in the far-N hotspot indicates a lower density, probably 
n e < 20 cm -3 . For the galaxy as a whole the integrated 
[SIIJ6717/6731 ratio is 1.191 (corresponding to 252 cm" 3 ). 

In NGC 4676, Wild et al. (2014) found intermediate 
[SII]6717/6731 ratios of 1.1-1.2 (n e ~ 200-400 cm" 3 ) in 
the bicones and both disks, without obvious radial gra¬ 
dients. The non-interacting M33 spiral, studied with the 
same PMAS spectrograph (Lopez-Hernandez et al. 2013), 
does show a gradient from an outlying HII region with 
[SIIJ6717/6731 ~ 1.5-1.8 (very low n e ), to the nuclear HII 


region with 1.30 (n e ~ 100), to the very centre with 1.07 
(n e — 328 cm" 3 ). In the starburst-merger Arp 299 (Heck¬ 
man et al. 1999) n e similarly increases from ~ 50 to 250-300 
cm" 3 at the centre. Le Tiran et al. (2011) found the stacked 
spectra of z = 1.3-2.6 starburst galaxies to indicate mean n e 
values of 200 cm" 3 (range 100-500) for the Ha-brightest re¬ 
gions, similar to local starbursts. On this basis, the n e ~ 750 
cm" 3 in the NGC 5394 nucleus is a fairly extreme case, prob¬ 
ably a consequence of the interaction. 

Kewley et al. (2013) predict in stellar photoionization 
models that increasing n e from 10 to 1000 cm" 3 could, by in¬ 
creasing collisional excitation, raise the N2 ratio by as much 
as 0.2 dex (their Fig 2; note that 03 is also increased by a 
similar amount). This might in itself account for the central 
high peak of N2 in NGC 5394, even without shocks or AGN. 
Observationally a high n e is associated with a positive N2 
offset in the BPT diagram for SDSS galaxies (Brinchmann, 
Pettini & Chariot 2008) and for individual HII regions in 
CALIFA (Fig 2 of Sanchez et al. 2015). Plotting for the in¬ 
dividual pixels in the 4 main star-forming regions (Fig 18), 
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NCG 5394 (spectra of two outer regions) 



Figure 19. Summed spectrum of the annular high-excitation re¬ 
gion where the N2 ratio exceeds 0.65, in total 296 arcsec 2 . Also 
shown for comparison is the (scaled-up) spectrum from a 14 pixel 
region of the S tidal tail. 


N2 is closely correlated with the [SII]6717/6731 ratio, which 
supports n e as the primary cause of the central enhance¬ 
ment. 

In contrast, the outer, annular region of high N2 and 
other ratios is not associated with any increase in n e . To 
confirm this region is genuine and investigate further we 
sum the spectra of all 296 pixels at 9 < r < 16 arcsec with 
N2 > 0.65 to give the spectrum shown in Fig 19. Ha and 
[Nil] 6584 are well-detected, with a ratio that appears close 
to unity on the plot, but in the integrated line fluxes from 
P0RT03D Ha is corrected upward for stellar absorption and 
the resulting mean N2 for these pixels (weighted by flux 
rather than area) is 0.741. In the same way the other mean 
ratios are found to be [OIII]/H/3 = 0.844, [SII]/Ha = 0.432 
and [OI]/Ha = 0.060. The [SIIJ6717/6731 ratio is 1.37, giv¬ 
ing a low n e of 20 cm -3 . The Ha emission and H<5 absorption 
EW are both quite moderate at 4.7A and 4.lA so this is nei¬ 
ther a starburst or a predominantly post-starburst region. 

The elevated line ratios in this region resemble (but 
are a little more moderate than) the shock-dominated outer 
regions of the late merger NGC 3256 (Rich et al. 2011), the 
extensive eastern shock region of the colliding NGC 7318 
in Stephan’s Quintet (Rodrfguez-Baras et al. 2014) and the 
bicones in the Mice. All 4 ratios would fit the models of Rich, 
Kewley & Dopita (2011) for composite shock excitation and 
HII regions, where the shock component is ~ 50% (of the 
Ha flux). Further comparison with the models tabulated by 
Allen et al. (2008) finds a dominant shock component that 
can fit the observed line ratios quite well with shock velocity 
v ~ 225-275 km s _1 and little or no pre-ionization (which 
would give a higher 03); this assumed near solar metallicity, 
magnetic equipartition (B = 3.23/xG) and a low electron 
density n e = 1 cm -3 . 

The total Ha flux (uncorrected for dust) from this shock 
region is 9.4 x 10~ 15 erg cm _2 s _1 (L = 3.2 x 10 39 erg s _1 ) 
only 2.6% that of the whole galaxy, so shocks remain a small 
(few percent) component of the galaxy’s emission-line lumi¬ 
nosity compared to star-formation. Note that this region 
only appears annular in ratio maps, while all the line fluxes 
increase towards the centre, so it is possible this emission 


component extends over the whole disk, but is not apparent 
inward of r ~ 4 kpc where star-formation dominates the Ha 
flux with a lower N2. In e.g., the elliptical NGC 4696, which 
has shocks from an ongoing merger but very little star for¬ 
mation (Farage et al. 2010), the N2 ratio remains uniformly 
close to 2. 

Similar high line ratios have also been observed in some 
passive, non-interacting galaxies, e.g. a mean N2 ratio 2.2 in 
the early-type CALIFA sample of Papaderos et al. (2013), 
which may be the result of hard continua from post-AGB 
(Asymptotic Giant Branch) stars, planetary nebulae or ‘hot 
low mass evolved stars’, but in these cases the lines are 
fainter, with EW < 2A. Fig 19 shows for comparison a 
scaled-up spectrum from a region in the southern tidal tail 
(also on Fig 4). The spectra are quite similar except the tail 
is a little bluer (less dust?), has stronger Balmer absorption, 
but lacks the emission lines (there is possibly a weak [Nil] 
line but the EW is no more than 0.8A). This is the part of 
the galaxy most dominated by the intermediate age stellar 
population. These stars cannot be making more than a very 
small contribution to the high excitation region, and we look 
elsewhere for the power source (see section 6.3). 


5 METALLICITY GRADIENTS AND THE BPT 

DIAGRAM 

Metallicities of galaxies or sub-regions can be estimated by a 
great variety of methods. One of the simplest requires only 
03 and N2 ratios to estimate the nebula gas metallicity, 
and can be applied even to high redshifts, with dust having 
little effect, but these ratios are subject to other influences 
(especially AGN) and there are uncertainties in calibration. 
The 03N2 index, given as logl0[O3/N2], is -0.387 at the 
centre of NGC 5394 and -0.184 for the galaxy as a whole. 
In the widely used PP04 (Pettini & Pagel 2004) calibration, 
metallicity expressed as the logarithmic ratio of oxygen to 
hydrogen, 12 + log(0/H) = 8.73 — 0.32 x 03N2. Marino et 
al. (2013), comparing line ratios from the CALIFA survey 
with calibrations based on electron temperature, derived a 
shallower relation, 12 + log(0/H) = 8.533 — 0.214 x 03N2. 

A different metallicity measure is obtained from the 
stellar component. Stellar and nebular metallicity may dif¬ 
fer, as stars preserve a record of the abundances at the time 
of their formation, which may be Gyr ago, and gas may 
move in or out. PORT03D model fitting generated maps of 
luminosity-weighted and mass-weighted stellar metallicity 
(lwZ and mwZ), the former giving more weight to recently 
formed stars, in units of Z (total mass fraction of N a tomic > 3 
elements), which are here converted to O/H on the basis that 
solar metallicity is Z = 0.015 and 12-|-log(0/H) = 8.72 (Al- 
lende Prieto, Lambert & Asplund 2001). Fig 20 shows the 
03N2 nebular metallicity averaged in radial annuli, using 
only pixels with F(Ha) > 10~ 16 and signal/noise > 3 for 
both ratios (so there is a gap between the disk and the N 
hotspot). Pixels with N2 > 0.7 or 03 > 0.7 are also excluded 
on the basis that they are likely to be more dominated by 
shocks than by star-formation, but this has little effect on 
the mean 03N2 values (shocks might cause O/H to be un¬ 
derestimated by increasing the 03 ratio but this may largely 
be compensated by the increase in N2). Fig 20 also shows the 
stellar lwZ and mwZ in the same annuli (using all pixels). 
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Figure 20. Radial profiles of the nebular (03N2) metallicities 
with PP04 and Marino et al. (2013) calibrations (with our fit¬ 
ted gradient), and of the luminosity-weighted and mass-weighted 
stellar metallicities from the PORT03D fits. 


Figure 21. ‘BPT’ diagram plot of the log of 03 vs. N2, 
for the pixels in the regions of strong star-formation, showing 
the local star-forming sequence from Kewley et al. (2013), the 
HH/composite demarcation of Kauffmann et al. (2003) and two 
iso-metallicity lines based on the Marino et al. (2013) calibration. 


The PP04 calibration gives a much higher metallicity 
than Marino et al. (2013), 12 + log(0/H) ~ 8.79/8.85 inte¬ 
grated/central (1.17-1.35Zq), compared to 8.57/8.62 (0.71- 
0.79Zq). The latter is much closer to the model-fit stellar 
mwZ. From r = 0.5 to r = 11.5 arcsec the Marino et 
al. 03N2 metallicity decreases by 0.080 ± 0.006 dex, but 
beyond this remains flat (perhaps even increasing slightly 
in the outer hotspot). A linear fit at r < 12 arcsec gives 
8.6166 —0.0069(±0.0002)r (plotted on Fig 20), and for r < 8 
arcsec only, 8.6205 — 0.0083(±0.0003)r. These gradients are 
-0.028 and -0.033 dex/kpc, or with the disk exponential 
scale h = 6.5 arcsec (Puerari et al 2005) corresponding to 
r e ff = 10.9 arcsec (2.75 kpc), can also be expressed as -0.076 
and -0.091 dex/r e //. 

Sanchez et al. (2014) found CALIFA disk galaxies typ¬ 
ically had a metallicity gradient of -0.1 dex/r e // at 0.3 < 
r < 2 r e ff, with flattening at even larger radii and evidence 
of a shallower mean gradient (-0.05) for interacting than for 
non-interacting (-0.11) subsamples (these estimates used the 
PP04 calibration and would be reduced by 1/3 for that of 
Marino et ah). Rich et al. (2012) found similar gradients, 
which flattened as mergers progressed. NGC 5394 appears 
consistent with their mean for early-stage mergers, but also 
with the non-interacting disks of the same mass in Ho et al. 
(2015). Simulations predict that interactions reduce the cen¬ 
tral O/H and flatten gradients by the infall of low-metallicity 
gas (Rupke, Kewley & Barnes 2010), but there is also a com¬ 
peting effect of re-enrichment by star-formation (Torrey et 
al. 2012), so the resulting evolution is more complex (in disk 
galaxies with high gas fractions the central O/H may even 
be increased). 

A new study by Barrera-Ballesteros et al. (2015), com¬ 
paring interacting and isolated disk galaxies in CALIFA, 
suggests that the effect of interactions on O/H is more 
complex than a simple progressive flattening. They found 
interactions had little effect on the average central O/H, 
and might even steepen metallicity gradients at small radii 
(r < 0.5 r e ff) where interacting galaxies show an enhance¬ 
ment in SFR. At larger radii (r ~ 1-2 r e ff) the mean metal¬ 
licity is reduced and it must be here the gradient is flattened. 


This was attributed to the inflow of lower metallicity gas. 
Also gradients in the outer galaxy could be ‘stretched’ out 
by the formation of tidal tails (Rich et al. 2012). These pro¬ 
cesses could explain NGC 5394 having a ‘normal’ O/H gra¬ 
dient at r < 10 arcsec which becomes flat only at r > r e ff- 

The mwZ shows no significant gradient in the disk (from 
r = 0 to r = 8.5 arcsec changing by —0.0018 ± 0.0331), 
whereas the lwZ falls steeply (—0.1135 ± 0.0151 between 
r = 0.5 and r = 6.5), and remains much lower than the 
mwZ in the outer disk, e.g. by 0.1359 ± 0.0222 at r — 8.5 (a 
6 a difference). This is surprising in view of the expectation 
that metallicity would increase with time, which would give 
lwZ > mwZ. Again a possible explanation is that the inter¬ 
mediate age and recent starbursts were fuelled by inflowing 
gas of lower metallicity compared to the underlying older 
disk stars, and, as Fig 8 of Barrera-Ballesteros et al. (2015) 
suggests, the resulting O/H reduction is most apparent not 
at the centre but at r ~ 5-10 arcsec (0.5-1 r e //). 

Returning to the nebular metallicity, Fig 21 shows the 
BPT (Baldwin, Phillips & Terlevich 1981) diagram, i.e. log 
03 vs. log N2, for the principal star-forming regions, to¬ 
gether with two isometallicity lines (using the Marino et 
al. calibration), and the ‘local star-forming sequence’ from 
Kewley et al. (2013), which is obtained from a photoion¬ 
ization model (with low electron density) and fits the lo¬ 
cus of typical z < 0.1 star-forming galaxies in the SDSS. 
The latter is essentially a sequence of increasing metallic¬ 
ity (right/downwards) and approximately orthogonal to the 
isometallicity lines [on which A(log N2) ~ A(log 03)]. The 
distribution of CALIFA pixels shows the increase in metal¬ 
licity from the outer star-forming regions to the central pixel, 
but combined with this, a large shift in the orthogonal direc¬ 
tion, along the isometallicity lines and away from the local 
sequence, by as much as 0.2 dex in N2. This latter cannot be 
due to metallicity but is as expected for a steep increase in 
Tie towards the centre (Kewley et al. 2013) and we propose 
this is the cause. The entire nuclear region transgresses the 
Kauffmann et al. (2003) demarcation, which was intended 
to divide pure star-forming galaxies from HII/AGN com- 
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Figure 22. Maps of the line-of-sight relative velocity A(u), from 
the red/blue-shifting of the emission lines Ha (left) and [NIIJ6584 
(pixels shown white, which have little emission, do not have valid 
measurements). Note they are closely similar. 


NGC5394: Emission-line Velocity Curve for PA=26 deg 



Figure 23. Line-of-sight velocity curve from the Ha wavelength 
averaged in ±15° sectors on the kinematically determined axis, 
position angle (f>k = 26° (East of North), with the NNE direction 
denoted positive. 

posites, but is too conservative in not fully allowing for the 
potential [Nil] enhancement of high n e and hard radiation 
fields, which can move HII regions across the line (e.g. Fig 
2 of Sanchez et al. 2015). 

6 KINEMATICS 

6.1 Rotation of the Disk 

Fig 22 shows a map of the line-of-sight velocity as measured 
from the Ha line, and also for [NIIJ6584. The galaxy is in¬ 
clined only 15° from face-on, according to K02, so that the 
line-of-sight A(u) is lower than the true rotation of the disk 
by a factor 3.9. Even so, the velocity gradient from rotation 
is clearly visible in the Ha velocity map. We obtain a velocity 
curve Vrot{r) by averaging this map in two opposite sectors 


Figure 24. Map of the stellar line-of-sight relative velocity A(i>) 
measured by fitting to the spectra (i.e. from absorption lines) on 
the V500 and the V1200 data. 


NGC5394: Stellar line-of-sight Velocity Curve for PA=26 deg 



Figure 25. Line-of-sight velocity curve from the stellar spectra, 
shown for both the V500 and V1200 data averaged in ±15° sectors 
on the position angle <Pt- = 26°. 


centred on the nucleus, in 2-arcsec bins of radius and within 
position angles 6 = 0fc±15 o and 6 = (</>*, +7r)± 15° where (j>k 
is the kinematic position angle of the galaxy. The velocity 
curve is evaluated in this way for all PAs <f> and rf)k esti¬ 
mated as the angle which maximises the velocity difference 
between the opposite sides, i.e. Sr=o 2 [ v ro t{r ) — v ro t{— r)]; 
this gave (j)k = 26±13°, in agreement with the kinematic (re¬ 
ceding side) PA of 25.8° measured for this galaxy by Garcia- 
Lorenzo et al. (2015). 

We take (j>k as 26° and Fig 23 shows the Ha v ro t (r) along 
this axis (NNE is positive, SSW is negative). The velocity 
curve shows the gradient produced by the galaxy’s rotation, 
with an amplitude (taken as half the difference between the 
mean of the two largest velocity measurements on each side) 
35 km s -1 . It is relatively symmetric, but with minima on 
both sides at r ~ 10 arcsec followed by steeper gradients, a 
feature seen in some barred galaxies and which might corre¬ 
spond to the weak distortion K99 found in the Ha velocity 
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map some 8 arcsec NE of the nucleus. The [Nil] velocity 
curve is very similar and gives an amplitude 43 km s -1 . The 
disk rotation velocity must be at least ~ 40 x 3.9 = 156 km 
s _1 , which is typical for a spiral galaxy in this stellar mass 
range (~ 2.4-4.6 x 10 10 Mg), e.g. Cortese et al. (2014). 

Shocks and ionized outflows are associated with mergers 
(Rich et al. 2011; Rich, Kewley & Dopita 2014), and might 
be revealed by kinematics, e.g. Yuan et al. (2012) attributed 
a high-N2 region to shocks/outflows as it was ~ 100 km 
s' 1 blueshifted with respect to the galaxy. The A(v) map 
for [NII]6584 might be more sensitive to shocks than Ha. 
However, the two velocity maps are very similar, and neither 
shows any region of blueshifted emission (ionized outflowing 
gas), either at the nucleus or in the outer (annular) high-N2 
region. 

Fig 24 shows the stellar A(v) maps, obtained by fit¬ 
ting the emission-subtracted spectra with STARLIGHT stel¬ 
lar models, for both the V500 data and the higher-resolution 
V1200, and Fig 25 the stellar velocity curves measured in 
the same way on the PA <f>k = 26°, centring on the nucleus 
in each case. The stellar A(v) maps, especially from V500, 
show a localized blueshift (~ 20 km s -1 ) near r ~ 0 and the 
cause is discussed below. The V500 stellar velocities show a 
weaker velocity gradient, with an amplitude (excluding the 
centre) only 21 km s^ 1 , but those from the V1200 data have 
a stronger gradient and a positive offset with maximum am¬ 
plitude 49 km s _1 . Cortese et al. (2014) found that stellar ro¬ 
tation velocities in SAMI (Sydney-Australian-Astronomical- 
Observatory Multi-object Integral-Field Spectrograph sur¬ 
vey) disk galaxies averaged 0.14 (±0.11) dex lower than v ro t 
from nebular emission lines, but considering the large un¬ 
certainties and low inclination it is not clear if the two differ 
significantly in this galaxy. 

6.2 Evidence of an Outflow 

We did not see evidence of outflowing ionized gas in the Ha 
and [Nil] velocity maps. Also both line profiles appear quite 
symmetric (although asymmetry was reported by K99), at 
least at the resolution of the V500 data, and are well-fitted 
by Gaussians. Any blueshifted line emission must be very 
faint compared with the star formation is the disk. How¬ 
ever, the cool gas phase of outflows may instead show itself 
in absorption lines, in particular MgII(2796,2803) and Nal 
(5890,5896). The nuclear spectrum plotted in Fig 4 has a 
very strong absorption line at 5890A restframe, which can 
be identified as the Nal (5890,5896) doublet, unresolved by 
the V500 grating and slightly blueshifted. While this line 
appears in purely stellar spectra (especially from K type 
stars), a stronger and blueshifted Nal absorption, seen in 
many LIRGs (e.g. Heckman et al. 2000), is a signature of 
starburst-driven outflows. 

To investigate this we sum the stellar model fitted by 
PORT03D over the same 5 pixels of the nucleus and compare 
this with the observed spectrum (Fig 26). The Nal line ob¬ 
served here with A = 5889.62 ±0.09A and EW = 7.03 ±0.06 
is > 3 times stronger than, and blueshifted with respect to, 
the purely stellar model where the blended doublet is fit¬ 
ted by a Gaussian with A = 5892.73 ± 0.24A with EW = 
2.26 ± 0.07A. The observed line blueshift is 3.11 ± 0.26A or 
158 ±13 km s" . Taking the (observed-model) residual spec¬ 
trum isolates the non-stellar, interstellar-gas component of 


the Nal line, which is here A = 5888.45 ± 0.09A, and by 
comparison with the (STARLIGHT- model) stellar Nal gives 
the blueshift of the non-stellar component as 4.38 ± 0.26A 
or 223 ± 13 km s -1 . 

This is indicative of a neutral gas outflow from the 
nucleus, directed perpendicular to the disk. In Fig 7, the 
spectrum of the nucleus region, integrated over 50 pix¬ 
els (r < 4 arcsec), shows enhanced Nal absorption with 
EW = 5.52±0.02A at A = 5890.12±0.02A (almost as strong 
as for the central 5 pixels only), but elsewhere the Nal line 
is consistent with the stellar models in EW (only 1 -2.5A) 
and A(~ 5892-3A). This confirms the strong blueshifted ab¬ 
sorption originates in the galaxy centre and is not a spurious 
feature produced by sky lines or some instrument problem. 
The outflow appears to have a similar, but weaker effect on 
the Call ‘K’ absorption line at 3933A; in the central 5 pixels 
the observed and stellar model wavelengths are measured as 
3931.48A (EW 3.6A), 3933.04A (EW 2.0A), and the resid¬ 
ual (excess absorption) is at 3930.30A, which corresponds to 
a very similar outflow blueshift A(v) = —209 km s -1 . This 
outflow might also explain the stellar velocity maps appear¬ 
ing slightly blueshifted at the centre. These A(u) could have 
been locally biased (up to ~ 20 km s -1 blueward) by the 
effect of the outflow on some absorption lines, especially in 
the V500 spectra which contain the Nal doublet. 

We calculate (integrating over wavelength) a map of 
the Nal EW (Fig 27); the regions of excess absorption ap¬ 
proximately fills the central r < 4 arcsec (r < 1 kpc) 
and, measuring with IRAF ellipse, has a FWHM 6.4-7.0 
arcsec (much larger than the instrumental FWHM of 3.5 
arcsec) and in the outer isophotes is slightly extended to 
the SE (at the EW = 3.5A isophote, PA = —36 ± 3 deg 
and ellipticity= 0.24 ± 0.02; using IRAF ellipse). Alongside 
we show a A(v) map obtained by Gaussian-htting the Nal 
absorption in each pixel (with IRAF htprof, over 5879 to 
5920A, avoiding the weak Hel emission line at 5876A). The 
blueshifted and high-EW regions correspond spatially, being 
extended by > 1.5 kpc on a PA close to the inner-disk minor 
axis, which seems consistent with an outflow perpendicular 
to the disk. Over this central region the total flux absorbed 
in the residual (non-stellar i.e. outflow) line is 1.44 x 1CP 14 
erg cnC 2 s _1 . 

The EW and A(v) of the Nal absorption are compa¬ 
rable to the Nal lines in many of the starbursting, Lir = 
101 ()_ i 2 Lq galaxies in Heckman et al. (2000), which have 
dust-rich super-winds. Blueshifted Nal absorption is ob¬ 
served to a lesser degree in more typical disk galaxies (Chen 
et al. 2010); the detection probability increases with central 
SFR, face-on orientation and high dust reddening, and the 
correlation of A(«) with disk orientation implies the outflows 
are perpendicular and entrained within bicones of opening 
angles < 60 degrees (Wild et al. 2014 measured this angle 
directly for NGC 4676A as n„,/47r = 0.37, corresponding to 
51°). The Nal blueshift/outflow velocity is also correlated 
with SFR, with ~ 200 km s _1 for 10 Mgyr -1 , and Bordoloi 
et al. (2014) found the same correlations for 2 ~ 1 galaxy 
outflows observed in Mgll 2796,2803. Our A(v) is slightly 
above average for the SFR but within the typical range. 

Heckman et al. (2000) give a formula (eq 5) for estimat¬ 
ing outflow rate as dM/dt ~ 60(r*/kpc)(Nu/3x 10 21 cm -2 )x 
(An/200 kms -1 ) x (fi ra /47r)M0yr _1 . Putting in starburst 
radius r* = 0.75 kpc, Nh = 6.1 x 10 20 at the centre (from 
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Figure 26. Observed spectrum of the central 5 pixels compared 
with the STARLIGHT stellar model fit, for (a) the whole spectrum 
and (b) for the region of the Nal(5890,5896) absorption doublet, 
together with the observed - model residual, showing that there is 
substantial blueshifted, non-stellar Nal absorption at the central 
region of NGC 5394. (c) This blueshifted absorption may also be 
visible, less strongly, in the Call (K) line at 3933A. 


Nal Absorption EW (Angstroms) Nal Absorption A(v) km s' 1 



pixels (arcsec) 


pixels (arcsec) 


Figure 27. (a) Map of the equivalent width of the Nal(5890,5896) 
absorption doublet, obtained by integrating over the line (5879— 
5904A), clearly showing the central region of strong absorption 
(up to EW = 7A), in excess of that expected for stars alone (about 
2A); (b) map of the velocity shift A(u) of the Nal absorption line 
relative to the stellar-fit 5892.73A, as obtained by a Gaussian fit 
to each pixel over 5879—5920A. 


Fig 5 of K02), Av = 223 km s -1 and S2 W /4 tt ~ 0.37 , gives 
dM/dt ~ 3.77 Mgyr -1 , similar to our estimate for the SFR. 
Heckman et al. (2000, and similarly Bordoloi et al. 2014) 
concluded that typically a ‘superwind is expelling matter 
at a rate comparable to the star formation rate’. The ki¬ 
netic power would be §^jjfv 2 = 5.9 x 10 40 erg s -1 ; 1/13 the 
galaxy’s dust-corrected Ha luminosity (Table 1). 


6.3 Evidence of Inflow? 

The outer high-excitation region (at 9 < r < 16 arcsec) does 
not show the same blueshifted A(v) as the outflow from the 
nucleus. Although very apparent in line ratio maps it seems 
to leave no imprint on the velocity maps in either Ha or 
[Nil], suggesting little motion perpendicular to the disk of 
the underlying galaxy. Nor is it associated with a signifi¬ 
cant broadening of the [Nil] line (at least at this resolution), 
which remains at FWHM ~ 6A across the galaxy. This may 
be similar to the late-stage merger NGC 3256, in which the 
outer regions with high, shock-like line ratios are blueshifted 
only ‘a few 10 km s -1 ’ (Rich et al. 2011), whereas the central 
strong Nal absorption (caused by outflow) is blueshifted by 
A(v) = —309 km s -1 (Heckman et al. 2000). 

The non-visibility of the outer high-excitation region in 
A(v) maps would argue that it is not outflowing gas driven 
by the nuclear starburst. We suggest that the starburst- 
driven outflow remains collimated in front of the central 
disk, as on the Nal map (Fig 27), and the outer high ex¬ 
citation region is powered instead by the kinetic energy of 
gas spiralling inwards within the disk plane, where it would 
have little A(v) in the line-of-sight direction, as a result of 
tidal forces from the interaction. The gas would lose angular 
momentum in the shocks and eventually form stars. 

This mechanism was argued for the outer regions of 
high excitation/line ratios found in a number of interact¬ 
ing or merging LIRGs (Monreal-Ibero et al. 2010, Rich et 
al. 2011). One argument was that the occurrence of outer 
high-excitation regions was more correlated with (late) in¬ 
teraction stages than with high SFR.. In the Torrey et al. 
(2012) merger simulations for a similar mass galaxy, the nu¬ 
clear gas inflow rate is generally ~ 0.5-1.0 x SFR and 1-5 
Mgyr -1 . At the velocities of disk rotation ~ 160-170 km 
s' 1 or as estimated for the shocks, ~ 250 km s -1 , and an 
inflow rate of at a few Mgyr -1 , the kinetic energy input 
would be sufficient to power the shock component of line 
emission (~ 10 4t> erg s -1 ). 


7 HISTORY OF THE INTERACTION 

NGC 5394 is now undergoing a nuclear starburst (generat¬ 
ing an outflow), has other localized knots of star-formation, 
and shows evidence of a more extensive (over most of the 
galaxy) episode of star-formation 10 s to 10 9 years ago. Fig 
28 shows the the PORT03D model-fits mass distributions of 
young (< 100 Myr) and intermediate age (0.1-1.0 Gyr) stars. 
The young population has a total mass 206 million Mg and 
is very centrally concentrated with 132 million Mg in the 
central < 4 arcsec and only 1.7 million in the far-N hotspot. 
The dust-corrected Ha SFR in Table 1, especially if fur¬ 
ther corrected to the Salpeter IMF of these masses (giving 
6 Mgyr - ) implies the current SFR is significantly higher 
than the mean (2.06 Mgyr -1 ) for the past 10 8 yr. Similarly 
our SFR of 3.39 Mg is more than twice the 1.54 Mgyr -1 
estimate of Lanz et al. (2013) based on SED-fitting with 
a model of the same IMF and 100 Myr of constant star- 
formation. Ha traces only the very recent SFR (< 10 Myr), 
compared to longer timescales for the infra-red emission, so 
the implication is again the SFR has increased strongly over 
the past 10 8 yr. 
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Figure 28. Distribution of the young (age zero to 100 Myr) and 
intermediate age stellar population (age 100 Myr to 1 Gyr), as 
estimated from the STARLIGHT model fits, shown on a log scale 
of stellar mass per pixel. 


The intermediate-age population has a very different 
distribution, in the outer disk as a broad ring and, especially, 
throughout the outer tidal arms. The earlier starburst must 
have been spread over the whole galaxy, whatever its shape 
at that time. The total mass of intermediate age stars is es¬ 
timated as 1.200 x 10 9 Mq with only 170 million M@ at r < 4 
arcsec. These patterns recall the merger simulation of Re- 
naud, Bournaud & Due (2015), in which early star-formation 
is extensive but the final burst is ‘almost exclusively central, 
due to gas inflows’. The concentration of star-formation in 
the SW disk (the inner western arm, see K02), reflects the 
lopsided distribution of Ht which would be a short-term ef¬ 
fect of the interaction, so is seen in the young population 
map but not for the intermediate-age stars, which instead 
are concentrated in the other two disk arms (eastern and 
outer-western). This agrees with to the ‘west-vs.-east time 
lag’ in star- and arm-formation (K99), with only the inner- 
western ‘third arm’ that is delayed i.e. recent. Similarly the 


two hotspots are not seen in the intermediate-age map, as 
they are younger features. 

K99/K02 describe NGC 5394 as ‘post-ocular’ and sim¬ 
ilar to the interacting galaxy IC 2163 but at a later stage 
of evolution, with star-formation more concentrated in the 
nucleus. IC 2163 is interacting with the larger spiral NGC 
2207 and indeed appears eye-shaped with most of its star- 
formation in numerous hotspots forming two arcs around 
the rim of the ‘eye’ (Kaufman et al. 2012). In the Struck et 
al. (2005) simulation of this pair there was a major starburst 
at the first close passage, some 280 Myr earlier, and there 
is predicted to be a second burst at least 100 Myr in the 
future as the galaxies finally merge. NGC 5394 could be 50- 
100 Myr ahead of IC 2163 and entering this stage, with the 
SFR increasing again as the galaxy approaches merger. In 
this simulation the long tidal arms are drawn out of the ini¬ 
tially round galaxy very soon after the first close encounter 
and starburst, and this might account for the tidal arms in 
NGC 5394 having the highest content of intermediate-age 
stars. 

The Arp 84 system is certainly at a later evolution¬ 
ary/merger stage than the Mice, for which the simulations 
of Wild et al. (2014) put the first perigalacticon only 170 
Myr ago, and observationally the oldest star clusters in both 
galaxies date from exactly that time (Chien et al. 2007). The 
Mice have less H<5 absorption and much less of an intermedi¬ 
ate age population than NGC 5394, and in the simulations 
of Wild et al. (2014) and Barnes (2004) their biggest star- 
burst lies 0.4-0.8 Gyr in the future when the galaxies will 
merge and the SFR peak at 4O-5OM0yr _1 . 

The hotspots in the northern arm of NGC 5394 are 
examples of the ‘hinge clumps’ seen in many interacting 
galaxies (Smith et al. 2014) - bright knots of active star- 
formation, typically found where tidal arms/tails join to the 
disk, with EW(Ho-) ~ 100A and L(Ha) ~ 10 4 ° erg s _1 . An¬ 
alytical and numerical models of tidal disturbance (Struck 
and Smith 2012, SS12) can explain many of the galaxy’s 
features including the formation of long, curved and paired 
tidal arms (which can extend to ~ 3 times the radius of the 
pre-interaction galaxy) and an inner disk with a structure 
described as ocular with internal arms (e.g. Figs 7 and 10 
of SS12); furthermore the model predicts ‘swallowtail caus¬ 
tics’ at the base of the tidal tails or arms where intersect¬ 
ing streams of gas would pile up, giving rise to localised 
star-formation which could be quite sustained and produce 
‘hinge clumps’ such as the inner-N hotspot. A refinement of 
the SS12 model, where a prograde encounter is represented 
by successive impulses, can produce additional swallowtail 
caustics far out from the disk on the spiral arms (Fig 12 of 
SS12) and in this way account for the formation and posi¬ 
tion of the far-North hotspot (and the similar feature in Arp 
82). Introducing a little asymmetry further improves the re¬ 
semblance to real Arp (1966) galaxies. The high-resolution 
simulation of Renaud et al. (2015) also predicts the forma¬ 
tion of a ‘young massive cluster’ far out on the northern tidal 
arm. Probably these hotspots will become globular clusters. 

In a further attempt to reconstruct the star-formation 
history we run STARLIGHT model fits on the 7 region spectra 
(Fig 7), with a set of Chabrier-IMF stellar templates repre¬ 
senting 23 ages from 1 Myr to 13 Gyr and four metallicities 
(Z=0.004, 0.008, 0.02, 0.05). For each f A (age, Z) template 
spectrum the program fits the mass of stars present now and 
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Figure 29. Reconstructed SFR history, estimated from the Star- 
burst model fits, of large (above) and and small (below) re¬ 
gions of NGC 5394 as a function of look-back time, plotted in 
A(log[age]) = 0.3 intervals. 

at formation. As the results are rather noisy, mean SFRs are 
estimated for wide A(log[age]) = 0.3 bins, by summing in 
each the (at-formation) masses of stars formed, and dividing 
by the bin width in years. In this model the current stellar 
mass is only 2.44 x 1O 1O M0, reflecting the ‘lightweight’ IMF. 
The model also fits dust extinction to the stellar continuum, 
A v = 1.06 mag for the nuclear region (notably less than the 
emission-line Au a = 2.32 mag from the Balmer decrement), 
and A v = 0.60 and 0.31 for the SW and outer disk. 

Fig 29 shows the resulting SFR histories. The nucleus 
has a large old (~ 10 Gyr) stellar population, plus resur¬ 
gent star-formation (note increasing SFR) over the past 100 
Myr. The outer disk and tidal arms are dominated by in¬ 
termediate age (10 9 -10 8 yr old) stars, the Ha-bright SW 
disk region has additional star-formation at age < 100 Myr, 
the star-formation in the far-N hotspot also goes back some 
100 Myr, while the inner-N hotspot, with stronger Balmer 
absorption, is mostly intermediate-age. The model fitted a 
large component of ~ 0.9 Gyr age stars for the outer disk 
and tidal arms, suggesting the interaction with NGC 5395 
and enhanced star-formation have been going on for a long 
time, perhaps even with two previous close passages. 

The first perigalacticon could have triggered a pro¬ 
longed influx of gas towards the nucleus, fuelling the sub¬ 
sequent star-formation. If this gas was relatively low metal- 
licity, the stellar metallicity would then decrease during the 
interaction, which (we argued in Section 5) could account for 
the (model-fit) luminosity-weighted stellar metallicity being 
lower than the mass-weighted (Fig 20). Yet this influx did 
not flatten the 03N2-metallicity gradient at r < 10 arc- 
sec, perhaps because of competing re-enrichment by star- 
formation (Torrey et al. 2012). Also 03N2 could be less sen¬ 
sitive to mixing of low-Z gas if prolonged star-formation and 
reprocessing produced a high central N/O ratio (e.g. Belfiore 
et al. 2015). Tidally-induced gas infall may be be the cause 
and the energy source of the relatively faint (a few percent 
of the total line emission) but extensive (~ 8 kpc diameter) 
outer ring of elevated N2, 03, [OI]/Ha and [SII]/Ha ratios, 
which are consistent with shocks (of velocity ~ 250 km s“ ). 
Extensive shock emission features similar to this this have 
been found in other late-stage or strongly interacting, gas- 


rich galaxies (Monreal-Ibero et al. 2010, Rich et al. 2011, 
2014; Rodrfguez-Baras et al. 2014). 

We also find evidence for a nuclear-starburst-driven out¬ 
flow of gas with (blueshift) velocity 223 ± 13 km s -1 , and 
estimated the rate as ~ 3.77 M 0 yr _ , which from Heckman 
et al. (2000) might not be exceptional if the SFR is near our 
estimate of 3.39 M 0 yr _1 . Wild et al. (2014) estimated 8- 
20 MgyW 1 at 350 km s -1 for the ionized outflow from NGC 
4676A. That the outflow in NGC 5394 is seen as absorption 
in the neutral line Na I, rather than in emission, might be 
a feature of the later evolutionary stage, with the neutral 
outflow persisting longer, as suggested by the detection of 
neutral outflows in some dusty, red post-starburst galaxies 
(e.g. Sato et al. 2009), and in recent gas-rich mergers (e.g. 
NGC 34 in Schweizer & Seitzer 2007). Much higher spectral 
resolutions (R > 5000) may be required to measure velocity 
dispersions and disentangle the different inflow, outflow or 
rotating components. 

Another feature of NGC 5394 is the high electron den¬ 
sity of n e ~ 750 cm -3 in the central kpc, compared to the 
central n e ~ 200-300 cm -8 more typical of star-forming 
galaxies. This is not unique, e.g. Krabbe et al. (2014), in a 
Gemini Multi-Object Spectrograph study of 12 galaxies in 
7 interacting pairs, found even lower [SII] 6717/6731 ratios 
and higher n e (~ 1000 cm -3 ) in some examples, (such as AM 
1054-325B and AN 1219-430), and concluded that there was 
evidence interacting galaxies on average had higher n e but, 
on the basis of moderate [OI]6300/Ha and [OIII]5007/H/3 
ratios, this was not generally caused by shock excitation. 
Otherwise no cause was given except to attribute a central 
increase in n e to gas inflow. The high central n e of NGC 
5394 could be a consequence of sustained gas infall over 
hundreds of Myr. Perhaps the n e enhancement is maximised 
for this particular interaction geometry (coplanar, prograde, 
with larger companion) and stage. Large integral field sur¬ 
veys which can measure line ratios diagnostic of n e will be 
able to quantify how it correlates with local and galaxy-scale 
properties. 


8 CONCLUSIONS 

(i) We examine CALIFA data for the spiral galaxy NGC 
5394, which is strongly interacting with the 4x more mas¬ 
sive spiral NGC 5395. Most of the Ha emission from NGC 
5394, and therefore the star-forming activity, is concentrated 
in the nucleus (75% in the central r < 1 kpc) with addi¬ 
tional star-forming regions in the SW of the disk and in two 
hotspots on the northern tidal arm. The galaxy is dusty, 
especially in the central r < 1 kpc where the Ha/H/3 ratio 
corresponds to a Ha extinction of 2.32 magnitudes. For the 
whole galaxy we estimate a dust-corrected, Ha based star- 
formation rate of 3.39 M 0 yr _1 (Chabrier IMF). The outer 
northern-arm hotspot has only 0.9% of the total SFR but 
the highest Ha equivalent width anywhere in the galaxy, 
over lOOA. Isolated, extranuclear concentrations of star- 
formation, like this hotspot, are predicted to form in merger 
simulations. 

(ii) Much of the galaxy outside the inner disk has 
a post-starburst spectrum with strong Balmer absorption 
(EW a bs(H<5) > 6A). The H<5 map and spectral model-fitting 
show the galaxy underwent an extensive episode of rapid 
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star-formation 0.2-1.0 Gyr ago. Comparison with simula¬ 
tions and other mergers suggest this was triggered by the 
earlier close passage of the two galaxies. 

(iii) NGC 5394 has a very high (for a purely star¬ 
forming galaxy) [Nil]/Ho ratio of 0.63 at its centre. The 
[OIII]/H/3, [OI]6300/Ho and [SII]/Ho ratios remain low 
here, with no central peak, and thus disfavour AGN and 
shocks as explanations. Instead the enhanced central N2 ra¬ 
tio is probably due to increased collisional excitation as¬ 
sociated with the very high central electron density, n e ~ 
750 cm -3 on the basis of the [SII]6717/6731 ratio. 

(iv) There is an outer region, forming an irregular ring 
of pixels around the nucleus at 9 < r < 16 arcsec (2.25 < 
r < 4 kpc), with even higher [NII]/Ha ratios of 0.65-0.85, on 
average 0.74. Here the [OIII]/H/3, [OI]6300/Ha and [SII]/Ha 
ratios are also high, suggesting a substantial (~ 50% of Ha) 
contribution from shocks (we estimate, with v ~ 250 kms -1 
and no pre-ionization). 

(v) The 03N2 index of nebular metallicity with the 
lower, Marino et al. (2013) calibration and the mass- 
weighted model-fit stellar metallicity are similar at 12 + 
log(0/H) ~ 8.6 (0.76 solar). The metallicity from the 
03N2 index shows a negative radial gradient of about -0.03 
dex/kpc (typical of disk galaxies) at r < 10 arcsec (or one 
rhi) , but this flattens at larger radii. Also the luminosity- 
weighted stellar metallicity is lower than the mass-weighted, 
especially in the outer disk, which might be the result of ear¬ 
lier inflow of lower metallicity gas. 

(vi) The Nal(5890,5896) absorption in the spectrum of 
the nucleus has an equivalent width 6 to 7A, at least 3 times 
stronger than expected from stellar absorption alone, and is 
blueshifted by 3A over a > 1.5 kpc region. The residuals 
between the observed line and a stellar model indicate the 
non-stellar absorption component is blueshifted by 223 km 
s . This is evidence for a neutral outflow of gas, at this 
velocity and perpendicular to the disk plane, driven by the 
central starburst. 

(vii) Although the galaxy is viewed only ~ 15° deg 
from face-on, a velocity gradient is visible in the Ha and 
[Nil] emission lines, and in the stellar spectra (at least for 
the V1200 data), consistent with a disk rotation velocity 
~ 156-180 km s _1 (as expected from the stellar mass). The 
Ha and [Nil] velocity maps are very similar and show no sig¬ 
nature of an ionized outflow (such as a blueshifted region) 
either at the nucleus or the outer ring of high N2 ratios. 
This suggests that shocks in the outer high-excitation region 
are not produced by the central outflow but more likely by 
interaction-triggered gas inflow within the disk plane. 

(viii) The strong neutral gas outflow, high electron den¬ 
sity and concentrated star formation ( > 1 MQyr -1 kpc -2 ) at 
the nucleus, and the strong Balmer absorption and hotspots 
elsewhere, are all a consequence of the prolonged interac¬ 
tion with NGC 5395. Over a timescale as long as ~ 1 Gyr 
and one or more perigalactica, this has reshaped the galaxy, 
formed tidal arms, triggered starbursts and gas inflow, and 
now the SFR may be increasing again as the two galaxies 
approach merger. 
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